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Strong mesonic decays are studied in the framework of the 't Hooft model for the 



two-dimensional QCD in the axial gauge. Special attention is payed to the processes 
with pions in the final state involved and the low energy theorems have been checked 
CnJ ' including the "Adler zero". It is demonstrated explicitly that in the chiral limit any 

three-meson decay amplitude with at least one pion vertex vanishes identically for any 
values of the external momenta. Ward identities for the dressed vector and axial- vector 
currents vertices are identified and their relationship to the pionic vertex are discussed. 

' 

q \ It has been recognized that strong decays play the crucial role in hadron spectroscopy 

issues. There are two reasons for this. First, as the distances involved in the light quark 
' sector are large, the dominant mechanism of strong decay should be nonperturbative, so that 
the decay studies deliver knowledge on the picture of confinement. Another point is that 
strong hadronic decays are proved to be a rather important tool in searches for constituent 
glue. The characteristic decay selection rules are known which distinguish between e.g. qq 
and hybrid mesons These selection rules are derived on the basis of quark models, and 
the question remains unanswered of how reliable are the model predictions from the QCD 
motivated point of view. 

In the constituent quark picture approach the confinement is modeled by a potential 
force, and the resulting spectrum describes the data reasonably well with notable exception, 
the pion. The chiral symmetry breaking (CSB) effects lie definitely beyond the scope of a 
simple constituent picture, and there is no hope to obtain Goldstone boson within any naive 
quark model. Similarly, decays with pions in the final state present problems for pair creation 
models. For example, in the standard 3 Pq model H, El |[] a reasonable fit to experimental 
data is obtained only with the so-called "mock meson" prescription for the phase space ||, 
while for the correct relativistic phase space the rate for processes with pions in the final state 
is too low [|J] . On the contrary, the microscopic model || of hadronic decays overestimates 
the decays with pions for higher quarkonia ||. What is worse, the naive pair creation models 
are based on the constituent picture and do not respect soft pions theorems; in particular, 
there is no hope to obtain Adler selfconsistency condition |7j for the amplitude with pions. 
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Unfortunately, the decays of immediate relevance for hadron spectroscopy are mostly the 
ones with pions in the final state ||. 

Adler selfconsistency condition follows from the most general symmetry considerations, 
and has nothing to do with the particular mechanism of confinement. Nevertheless, the 
idea that confinement and CSB phenomena are interrelated seems to be rather meaningful. 
A model was suggested many years ago || which connects confinement and CSB (see also 



e.g. IJTOj, where similar ideas were employed). Unfortunately, this model main ingredient, 
three-dimensional oscillator confining force, is not motivated by QCD, and among other 
drawbacks of such an approach is lack of gauge and Lorentz invariance. An important 



step was made in [HJ, where it was shown that the same gauge and Lorentz invariant 



nonperturbative gluonic correlators which produce area law are responsible for the formation 
of chiral quark condensate. Thus the quark model motivated by such a picture should be 
free of above-mentioned drawbacks and should be able to reproduce, inter alia, all pion 
properties. Obviously this quark model is going to be rather unusual, and it is instructive 
to study strong decays in some toy exactly solvable theory with confinement and CSB. 

An example exists of a theory where both confinement and CSB are due to the same force, 
and it is two-dimensional quantum chromodynamics (QCD2) in the large Nc limit which was 
first considered many years ago |l2j and still remains popular in studies of various aspects 
of strong interactions. The large Nq limit is important in establishing the chiral properties 
of the theory. Indeed, it is well known that the Coleman theorem |L3) prohibits CSB for 
any finite number of colours in a two-dimensional theory. Still there is no contradiction with 



the Coleman theorem ]T3|] if one considers the weak coupling regime of the theory where 
m q ^ 9 ~ l/V^C; 2.e. the limit of infinite number of colours is taken first (see e.g. |14j for 
detailed discussion of this issue as well as of another phase of the theory which corresponds 
to the strong coupling regime m q <C g and where the Coleman theorem works at full scale). 

Most of the studies in QCD 2 were performed in the light-cone gauge, which considerably 
simplifies the spectrum calculations but yields perturbative vacuum. There exists an alter- 
native approach based on the Coulomb gauge A\ — It appears that the vacuum is 
nontrivial in the Coulomb gauge version, and nonzero quark condensate exists for massless 
quarks [IB]. The latter feature is confirmed by sum rules approach |TJ| in the light-cone 



gauge. 

In the present letter we perform studies of strong hadronic decays in the QCD 2 employing 
the Coulomb gauge. This choice enables treating pions on the same footing as other mesons. 
The latter is in contrast to light-cone gauge choice, where the pion wave function is very 
singular, and one is forced to use sophisticated nonperturbative methods, like sum rules and 



operator product expansion, to arrive at reliable results in the pion physics (Tj| . 

Let us review the essentials of QCD 2 in the Coulomb gauge (see |15|] for the details). 
Our convention for 7 matrices is 70 = 03, 71 = ia 2 , 75 = 7o7i- The large Nc limit implies 
that g 2 Nc remains finite. The gluonic propagator in the Coulomb gauge A± — takes 
the form D 0Q (k ,k) = —1/k 2 , and the infrared singularity is regularized by the principal 
value prescription yielding the linear confinement. There are no transverse gluons in two 
dimensions, and in the large Nc limit only planar graphs survive, so the only nontrivial 
one-particle quantity is the quark Green's function 

S(po,p) = 1 w r , (1) 
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which can be found from the Schwinger-Dyson equation with the result 



E(p) = [E(p) cos 9(p) — m] + ji[E(p) sin 9(p) — p], (2) 

pcos8(p) — m sin 6(p) = — j- — sin[0(p) — 0(k)], (3) 

Zi J [J) rC J 

E(p) = m cos 6(p) +p sin 6(p) + - j- —cos[6(p) - 6(h)], (4) 



where 7 = 9 ^ c . Parameter 6(p) has the meaning of the Bogoliubov-Valatin angle describing 
the rotation from bare to dressed quarks. The obvious properties of the solutions to the 
equations (§), (@) are 6(-p) = -6(p), E{-p) = E(p), and 8(p) -^ir/2. 

The spectrum of the theory is defined from the homogeneous Bethe-Salpeter equation 
for the matrix wave function 

S(p, P) = T{p) (^p 7 5^ + (p, P) + ^^75^-(P, P)) T + (P - p) (5) 
with T{p) = exp(— ~0(p)7i). The functions <p±(p, P) are solutions to the system of equations 
[E(p)+E{P-p)-P Q ]ip + (p,P) 

/elk 
-—^ [C(p, k, P)tp + (k, P) - S(p, k, P)cp-(k, P)] 

(6) 

[E(p)+E(P-p) + P ]<p-(p,P) 

/A V* 
— [dp, k, p)<p-(k, p) - s{p, k, p)<p + % p)\ , 

where 

C{v , *, P) = cos cos 'C-ri -»(*-*>, 

S<p,k,n = *W-W*n« P -''>- 9 < P - k \ (7) 

with P and p being the total meson momentum and the momentum of the quark respectively. 

The system (§) describes a meson moving forward in time as a superposition of the 
quark-antiquark pair moving forward in time with the amplitude <p + and backward in time 
with the amplitude <£>_. Such a particle- hole interpretation is supported by the Hamiltonian 
approach developed in |T7[, where the operator m+(P) creating the n-th meson with total 
momentum P was constructed as a linear combination of operators creating and annihilating 
the qq pair. 

A very interesting point concerning Bars-Green system of equations (0) is that if it is 
treated as a matrix integral equation, then its kernel is not Hermitian fl7| . The eigenvalues 
Pq are real, but the orthonormality and completeness conditions take the form (see [I7| for 
details) 

J^(VI(P,P)V+(P,P)-V-(P,P)<P-(P,P)) = 

(8) 

J^(tpl( P ,P)^ip,P)-^_ip,P)^(p,P)) = 
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oo 

E(v+(p.^(*.^)-v-(p.^)^(*.^)) = o- 



(9) 



n=0 



Solutions of the system (^j) come in pairs: for each eigenvalue Pq with eigenfunction 
V 9 -) there exists another eigenvalue —Pq with eigenfunction <£>'"). With this sym- 
metry only positive eigenvalues enter the forms (§), ([|). 

The system (||) is solved numerically in and it is shown that the component is 
small i) if the quark mass is large and ii) for higher excited states. In both cases quark 
potential model with local linear confinement serves as a good approximation. Besides that 
the component dies out with the increase of the total mesonic momentum P: in the 
infinite momentum frame (P — > oo) it is zero and the equation for the p + reproduces the 
't Hooft equation [I2| after an appropriate rescaling, as it was shown in ]T5|] . 

The chiral condensate is calculated straightforwardly from equation ([I]): 



(vac\qq\vac) 



N, 



c 



2tt 



dk cos 9{k). 



(10) 



It is shown in [[L(| that the gap equation (Q) has nontrivial solution in the chiral limit 
m — 0, and the chiral condensate (10) does not vanish with this solution. It means that the 



Goldstone mode exists in the spectrum, and its explicit form was found in [17 



<Pl(P,P) = \l2p 



71 



cos • 



9(P-p)-0(p) 9(P-p) + 9(pY 



11) 



for P > and Pq = VP 2 ((fl <-+ <p\ for P < 0). For P = <^+(p, 0) = ip!(p, 0), pion spends 
half of time in backward motion of the pair, and such a function has zero norm, as it should 
be for the massless particle at rest. In the opposite limiting case P — » oo the backward 
motion part dies out, and 




< p < P, 



(12) 



coinciding with the Goldstone mode of the 't Hooft equation. Form (|12| ) is the main reason 
to consider the Coulomb gauge version of the theory, as all nontrivial content of the wave 
function (12) is concentrated in the boundary regions x — > and x — ► 1 (x = p/P). The 
same is true of course for the QCD 2 quantized at the light-cone [Bl. Quantities like chiral 
condensate do not come out trivially with such a singular wave function fHf . 

Now we are in the position to calculate the strong decay amplitude (Fig.l), which is 
0(l/y/Nc)- The answer reads 

M(A -> B + C) = 



2 Sp(T A (k + P B ,P A )S(k-P c )T c (k,P c )S(k)T B (k + P B ,P B )S(k + P B )) (13) 




4 




where the meson-quark-antiquark vertex Ym can be expressed in terms of the Bethe-Salpeter 
wave function (^) 

fdk $ M (k,P) 
rM(p,P) = J ^70^^70, (14) 

and r M (p, P) = 7 r| f (p, P)7 . As the vertex ([HI) does not depend on p > the integration 
along the loop in (|T^) is trivial, so one arrives at the following final expression for the decay 
amplitude in the initial meson rest frame Pa = 0, Pg = —Pq = P' 



M{A^B + C)= (15) 
{ — (ft(k + p, 0)(pH(k + p, 0) [c{-p, q, k)(p%(q, -p) + s(-p, q, k)(p- C ] 



Vn^J (q-k) 2 

-tfi (k + p, 0)<f%(k, -p) [c(p, q + p,k+ p)(p+{q + p,p) + s(p,q + p,k + p)(p^(q + p, p)] 
-tp%(k, -p)<p^(k + p, p) [s(0, q + p,k + p)(pA(q + p,0) + c(0, q + p, k + p)ip*{q + p, 0)] 
+ip c (k, -p)ifl{k + p, p) [c(0, q + p,k+ p)ifi(q + p,0) + s(0, q + p,k+ p)^(q + p, 0)] 
+(f A (k + p, 0)(p^(k, —p)[s(p, q + p,k + p)(p+(q + p,p) + c(p, q + p,k + p)(f B .(q +p,p)] 
+^>+{k + p, 0)<p+(k + p, p) [s(-p, q, k)(p+(q, -p) + c(-p, q, k)<p°(q, -p)] } 

+(B <-> C,p <-> —p), 

where 

c(p , q , k) = C0B m^M sin g(P-*)-»(P-ri 

Amplitude ( ^5|) could be easily obtained from the Hamiltonian approach developed in 
[|l7j], though further calculations appear more transparent using the matrix wave function 
(|5|) and the matrix form of the bound state equation (|6[). Nevertheless, equivalence of both 
approaches can be established at any intermediate step. 

If one neglects the backward motion contributions in the amplitude ( |To] ) and inserts the 
nonrelativistic values of the angle 9 (cos8(k) = 1, sin8(k) = k/m) then it reproduces the 
standard quark model decay amplitude due to OGE Coulomb interaction || adapted to the 
two-dimensional case. It is clear, however, that substitution of the nonrelativistic angle is 
not justified for kinematical reasons. 

The six-term form ( |T5P was anticipated in |[20|| , where the role of backward motion in 
decay was discussed for the first time. It was argued there that a simple procedure could 
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Figure 2: The meson- vector (a) and meson-axial- vector (b) current couplings. 



be used as a crude approximation. Namely, one neglects the backward motion pieces for 
all mesons but pions. Then a simple graph counting would give the net result: amplitudes 
with a pion in the final state are multiplied by an overall factor 2 over the naive quark 
model result, and those with two pions are multiplied by 3. Obviously this prescription of 
constructive interference does not solve the Adler selfconsistency problem. 

It appears that in the two-dimensional amplitude (p~5|) the interference between the six 
terms is destructive: the total amplitude vanishes in the chiral limit, if at least one 
of the final state mesons is the pion. Of course, one can perform explicit calculations 
inserting the pion w.f. ([□]) and making use of the equations (Q), (^) and (Q). There exists, 
however, a simple and elegant expression for the pion-quark-antiquark vertex in the limit 
m = 0: 



T n (p,P) = S- 1 (p)(l 



75 ; 



:i-i 5 )s-\ p -p). 



(16) 



Inserting it into equation (13) trivially yields M(A — > n + C) = in the chiral limit. 

Expression flTBD for the pion vertex follows from rather general considerations and is 
connected with the vector and axial-vector currents conservation. 

The meson-vector current coupling V?f{P) = {vac^^M, P) can be calculated as 
(Fig.2a) 



V. 



M 



ijjNc 



I T^ Sp{s{p ~ p M P )r M (p, p)) = JT c j ^.sp(j^ M (p, p)). (17) 



dp 



(2vr 



Substituting expression (151), one finds 



M 



-^J N c I ^[<P+ {p, P) + <P-(p, P)) sin 



6(P -p) + 6{p) 



V 



M 



To demonstrate the vector current conservation we write down the equation for the matrix 
wave function (|5T) as 



P M $mG»,P) = (KP+7om)$ M (p,P)-®M(p,P)(l5(P-p)+7om) 



+ 7 



dk 



(p — k)'- 



{A + (k)$ M (p, P)A_(P - k) - A+(p)$ M (fc, P)A_(P - p) 



(19) 



-A_(A;)$ M (p, P)A + (P - k) + A_(p)$ A f(A;, P)A + (P - p)} 



where projectors are defined as A±(p) = T(p)^T+(p). 
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Then, integrating this equation over p, multiplying by 70 and taking trace one gets 

pMyM _ pyM = 0> (2Q) 

i.e. the vector current is conserved. 

Similarly, for the axial- vector current matrix element A^f(P) = (vac\ 97^759! M, P) one 
can find (Fig. 2b) 

PfA? - PA M = -2my^ f pSp($), (21) 

that yields axial- vector current conservation for m = 0. Slightly relaxing the chiral limit one 
gets the pionic decay constant f n = J ^- and the Gell-Mann-Oaks-Renner relation [JH 

flMl = -2m(qq) (22) 

(see also |IT||). It is instructive to mention that in the chiral limit vector current has 
nonzero coupling only to the pion: the vector current matrix elements (|18|) contain the meson 
wave functions integrated with the pionic one (|TTJ) , so the corresponding decay constant 
vanishes due to orthogonality conditions (§). (We remind here that the pion does couple to 
the vector current in 1+1!). As the axial- vector current in 1+1 is dual to the vector current, 
the axial decay constants Jm for excited mesons vanish in the chiral limit too. 

Calculations of mesonic form factors (Fig. 3) are more involved, and require knowledge of 
the full quark- ant iquark scattering amplitude T ik i m (p^, k^, P M ), 

27rz^y ^ 1 — 

r ik ,im(Pn> K> p v) = 7 TT^{lo)ii{lo)k m - h 2 -5 5Ffjtf(Pi P)) im (T M (k, P) kh 

VP ~ K ) M ~~ ^0 

+ *7 2 E P ! p aA ^ T m(P ~ v)lo)imW M {P ~ k, P)j ) k i- (23) 

which is the solution of the inhomogeneous Bethe-Salpeter equation. 

The dressed vector current-quark-antiquark vertex (Fig. 4) is then given by 

1 f dq — 

v^p,Pv) = «7 M + «7Ep — m / 7r s Pb^M(qiP)]?M(j>,P) 

u -Ml _ -in J 



M 



'7E P 1pm I ^-Spij^UP-Q^^loTMiP-p^ho, (24) 
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Figure 4: Graphical representation of the dressed vector current-quark-antiquark vertex. 



where P is the total momentum of the quark- ant iquark pair, and p is the momentum of the 
quark. One easily recognizes the vector current couplings (|17D entering (|24|), so that in the 
chiral limit only pion contributes into v^, displaying a sort of "vector dominance" . The same 
takes place with the axial-vector current-quark-antiquark vertex. 

Using the equations ([!]), (fj) and (§) together with the completeness condition (H) allows 
to recast the vector current divergence into the Ward identity form (similar expression in 



the light-cone gauge was derived in [22 



- ii> M (p, P) = S-\p) - s-\ P - P), 

which makes vector current conservation manifest for the form factors: 

Q,(M, P\v„\M', P') = 0, Q, = P»- P;. 
Similarly, the axial-vector current divergence in the chiral limit is 

- iP^(p, P) = s-\ P ) l5 + l5 s-\ P - P), 

and 



(25) 
(26) 
(27) 

Q,(M,P\ afl \M',P') = 0, g M = P^-P;. (28) 

The pion-quark-antiquark vertex (|i~6|) is nothing but a linear combination of the vector 
and the axial-vector currents divergences (|25|) and ( p?|) 

r ff (p, P) = -iP^(p, P) - zi>„(p, P). (29) 

Relation (|29|) plays the role of PCAC hypothesis in 1+1: as pion couples both to vector 
and axial- vector currents, the form (|29|) is inevitable in 1+1 then. 

In conclusion, let us briefly recall the results reported in the present letter. In the previous 
publication |17| we developed a Hamiltonian approach to the two-dimensional QCD in the 
axial gauge and demonstrated how the Hamiltonian of the model could be diagonalized 
in the mesonic sector. One of the key features of the above Hamiltonian is the distorted 
norm (^) plugged into it which allows description of the massless qq states, chiral pions, on 
equal footing with the rest of higher massive excitations. Chiral properties of the model 
were reproduced and discussed. In the meantime a crucial role of the backward in time 
motion of the qq pair inside meson, extremely important for the pion, was anticipated but 
not discussed in detail. In the present letter we return to this issue and investigate the role 
played by pions in the strong hadronic decays and in providing the low energy theorems 
like "Adler zero" selfconsistence condition. We use an effective diagrammatic techniques 
involving dressed meson-quark-antiquark vertex, quark and mesonic propagators and the full 
four-quark scattering amplitude. Any hadronic process can be formulated diagrammatically 
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using the above ingredients. The three-meson amplitude is studied in detail in the suggested 
approach and a strict result is obtained stating that the interference between different parts of 
this amplitude which are due to (p + and y?_ components of the pion is completely destructive, 
so that the total amplitude vanishes for any pion momentum, not only in the limit p n — > 0. 
The latter result seems quite natural if one takes into account that the pionic decay constant 
f w is dimensionless in 1+1 and does not define any scale for soft pions then. 

Developing similar approach for the four- dimensional QCD, which could shed light on 
many puzzles of hadronic spectroscopy, is still an open challenge. In the model || the 
interaction is described by the time component of a vector force, yielding the gap equation 
similar to (|3|), (^) and the Bethe-Salpeter equation similar to ([19]) • The existence of a chiral- 
noninvariant solution of the gap equation implies the existence of a Goldstone boson. Axial- 
vector current is conserved in the chiral limit of the model, and all the relations of current 
algebra are satisfied. However, the model is not covariant, that prevents from proceeding 
further along the lines presented here. We point out once more that the approach suggested 
in [|ll|] is rather promising: the gauge invariant confining interaction employed in [|TTJ is 
Lorentz covariant by construction, respects the area law and is responsible for the chiral 
condensate formation. We expect that the quark model derived within such a formalism 
could be able to reproduce all the nice features described above and, being fully covariant, 
could allow to calculate hadronic observables. 

The authors would like to thank B.L.Ioffe for enlighting discussions. Financial support 
of RFFI grants 00-02-17836 and 00-15-96786 and INTAS-RFFI grant IR-97-232 is gratefully 
acknowledged. 
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